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Abstract 

Influence of the weak electric field on the electronic 
structure of the Fibonacci superlattice is consid- 
ered. The electric field produces a nonlinear dy- 
namics of the energy spectrum of the aperiodic su- 
perlattice. Mechanism of the nonlinearity is ex- 
plained in terms of energy levels anticrossings. The 
multifractal formalism is applied to investigate the 
effect of weak electric field on the statistical prop- 
erties of electronic eigenfunctions. It is shown that 
the applied electric field does not remove the mul- 
tifractal character of the electronic eigenfunctions, 
and that the singularity spectrum remains non- 
parabolic, however with a modified shape. Changes 
of the distances between energy levels of neighbour- 
ing eigenstates lead to the changes of the inverse 
participation ratio of the corresponding eigenfunc- 
tions in the weak electric field. It is demonstrated, 
that the local minima of the inverse participation 
ratio in the vicinity of the anticrossings correspond 
to discontinuity of the first derivative of the dif- 
ference between marginal values of the singularity 
strength. Analysis of the generalized dimension as 
a function of the electric field shows that the elec- 
tric field correlates spatial fluctuations of the neigh- 
bouring electronic eigenfunction amplitudes in the 
vicinity of anticrossings, and the nonlinear charac- 
ter of the scaling exponent confirms multifractality 
of the corresponding electronic eigenfunctions. 
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1 Introduction 

Statistical properties of the electronic states in 
nanosystems are a subject of great interest (see, 
for example, Refs. [TJ [5J [31 0] and the references 
therein) . Partially it is due to the progress of exper- 
imental methods of nanophysics which allows to in- 
tentionally fabricate high-quality heterostructures 
consisting of alternating layers of different materi- 
als [51 [71 [S] . The thickness of each layer can be 
controlled during the growth process with accuracy 
of one atomic monolayer, so that one can fabricate 
multilayer systems (superlattices) with the desired 
geometrical parameters of layers and well defined 
interfaces. In this way periodic as well as disor- 
dered multilayer systems can be obtained by the 
sequential deposition of layers with different thick- 
ness of material. 

An intermediate case between periodic and disor- 
dered multilayer systems corresponds to aperiodic 
order of layers [HI [TO] . Such structures are inten- 
tionally generated by the deposition of layers of 
two different materials according to the Fibonacci, 
Rudin-Shapiro, Thue-Morse, etc. sequences [TTTfT^l 
H3l HH [151 HH HH [HI HH EDI EU E21 ESI. Exper- 
imental as well as theoretical studies of these su- 
perlattices are concentrated on the consequences 
of the long-range correlations induced by the ape- 
riodic arrangement at a length scale longer than 
atomic one [HIES]. In particular, this problem has 
been extensively investigated in the Fibonacci su- 
perlattices which are regarded as a typical example 
of aperiodic systems [551 (2Z1 HE]- ^ n t nese stud- 
ies, it has been found that the wave functions of 
one-particle states are critical, i.e. nor extended 
neither localized [29]. Further studies have shown 
that the decay of the envelope wave function obeys 
the power law and its structure can be regarded 
as a multifractal resembling the case of electronic 



wave functions in disordered systems at the mobil- 
ity edge [I] [3U] . The shape of the wave function 
is highly fragmented in the finite Fibonacci super- 
lattice, and in the limit of infinite Fibonacci super- 
lattice corresponds to a self-similar Cantor set with 
zero Lebesgue measure [3Tj [35] . 

The purpose of this work is to examine the ef- 
fect of the nonlinear dynamics of the energy levels 
driven by weak electric field on the global as well 
as local electronic structure of the finite semicon- 
ductor Fibonacci superlattice made of two different 
semiconductor layers. Energy spectrum of this su- 
perlattice is particularly interesting in weak electric 
fields where anticrossings lead to closing of the en- 
ergy gaps [33J GUI ES] ■ Particular attention is paid 
to the effect of this nonlinear dynamics of energy 
levels on spatial fluctuations of the electronic wave 
function amplitudes and their spatial extents. 

In the systems with broken translational sym- 
metry the multifractal analysis provides deep in- 
sight into the nature of the electronic wave func- 
tion [311 133 EH1 EH1 HOI SH S21 H3] and allows to 
explore the localization phenomena in the presence 
of the electric field. 

Although the results of the present work are re- 
lated to the semiconductor superlattice, they can 
be also generalized to other aperiodic superlattices, 
e.g. photonic or phononic band gap structures un- 
der influence of appropriate perturbations [44 [ 145] . 

The paper is organized as follows. In Sec. [2] we 
present the model of the Fibonacci superlattice and 
the methods of its analysis, in Sec. [3] we present 
the results and discussion, and the conclusions are 
presented in Sec. [4] 

2 Model of the Fibonacci su- 
perlattice and the methods 
of analysis 

We consider the semiconductor aperiodic super- 
lattice generated according to the Fibonacci se- 
quence of two different semiconductor layers made 
of Alo.3Gao.7As and GaAs The differences 

in the bandgaps of Alo.3Gao.7As and GaAs semi- 
conductor layers lead to the discontinuities in the 
conduction as well as valence band edge profiles at 
the interfaces [15] . This creates the effective poten- 
tial energy that consists of the set of barriers and 
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Figure 1: (Colour online) Potential energy profile 
of the Fibonacci superlattice formed by 20 layers of 
Alo.3Gao.7As and 31 layers of GaAs for without the 
external electric field. The subsequent layers are 
composed following the Fibonacci binary sequence, 
i.e. each Alo.3Gao.7As or GaAs layer corresponds to 
1 or in the sequence, respectively. The electronic 
states are plotted as (red) horizontal lines at their 
energies, with lines plotted for those x at which 
exceeds its average value. 

potential wells distributed along the growth axis of 
the semiconductor superlattice (Fig. [T]) . 

The conduction-band potential energy V(x) is 
modelled by the superposition of the power- 
exponential potentials in the form |47j 
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K(x) = £V exp 

i=l 
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where N is the number of Alo.3Gao.7As barri- 
ers located at positions Xi and having height Vq, 
with parameters c and p characterizing the shape 
of barriers. The positions of barriers, x i , are 
distributed according to the binary Fibonacci se- 
quence generated over set {0,1} using the follow- 
ing inflation rules gS] EH]: — ► 01, and 1 — > 
0. These rules allow us to obtain the sequence 
(0, 1, 0, 0, 1, 0, 1, 0, 0, 1,0,.. .) of any desired length. 
In our notation 1 corresponds to Alo.3Gao.7As layer 
(barrier) and corresponds to a single GaAs layer 
having the same width as the barrier. 

The envelope wave function of one-particle elec- 
tronic eigenstate of the Fibonacci superlattice can 
be found by the solution of the Schrodinger equa- 
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Figure 2: (Colour online) The motion of the eight 
lowest eigenvalues (n = 0, 1, .., 7) driven by the elec- 
tric field in the finite Fibonacci superlattice formed 
by 20 layers of Alo.3Gao.7As and 31 layers of GaAs. 

tion within the effective mass approximation 

+ = ( 2 ) 

2m* dx z 

where m* is the effective mass of the conduction 
band electron. 

The Schrodinger equation ([2| with the potential 
energy V(x) given by equation and the Dirich- 
let boundary conditions, ip(0) = tp(L) — (where 
L is the size of the computational box), constitutes 
the eigenvalue problem which can be solved numer- 
ically. This form of the boundary conditions allows 
us to neglect the complexity of the energy spectrum 
due to the surface states. 

The energy spectrum of the Fibonacci superlat- 
tice formed from a finite number of barriers and 
quantum wells is characterized by the density of 
states which can be expressed as follows 

BOS(E) =Y,KE-E n ). (3) 

n 

This quantity gives only global properties of the 
system, whereas the local properties can be de- 
scribed by the local density of states which directly 
explores the local amplitude of the electronic wave 
function for a given energy. The local density of 
states (LDOS) is defined by the formula 

LDOS(x, E) = Y, \Mx)\ 2 5(E - E n ), (4) 



where ip n (x) is the normalized electronic wave func- 
tion of the n-th eigenstate. 

One of the most characteristic properties of the 
electronic states in the finite Fibonacci superlattice 
are the self-similar spatial fluctuations of the elec- 
tronic wave function amplitude which stem from 
the aperiodic distribution of the barriers and wells 
in the system. These fluctuations can be analyzed 
by the multifractal formalism which allows one a 
deeper insight into their nature. An essential in- 
gredient of the multifractal formalism is the nor- 
malized probability measure of the wave function 
in the fc-th box Bk of linear size e, 

P nk (e)= [ dx\^ n (x)\ 2 - (5) 

JB h {e) 

The multifractal analysis of the electronic wave 
functions in the Fibonacci superlattice can be per- 
formed efficiently by applying the box-counting 
procedure [5U1 151] with condition: a <C e <C L (a is 
the lattice constant). The normalized g-th moment 
of the probability measure of the wave function is 
given by the formula 

where N(e) = L/e is the number of boxes. 

For each value of the scaling index q we can de- 
termine the singularity strength according to the 
formula 

a n (q) = lim S^W^/We; 1) (?) 

and the corresponding singularity spectrum in a 
parametric representation as follows 

/K) - Inn — , (8) 

where S = e/L is the ratio of the box size to the 
system size. 

The singularity spectrum f(a n ) gives an accu- 
rate description of the multifractal properties of the 
probability measure of the wave function. Instead 
of the singularity spectrum, we can equivalently 
consider a hierarchy of generalized dimensions of 
the wave function. It stems from the fact that the 
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generalized dimension D n {q) is related to the sin- 
gularity spectrum by the formula 



D n (q) 



(9) 



For the integer values of the scaling index q, the 
generalized dimensions of the wave function have 
a physical meaning |52| . A particularly interest- 
ing generalized dimension corresponds to q = 2 
when it is known as the correlation dimension of 
the wave function. Its significance results from the 
relation to the density-density correlation function 
and the inverse participation ratio (IPR) which is a 
measure of the spatial extent of the eigenstate and 
can be used to describe the localization properties 
of the electronic eigenstates in real space [S3]- It 
is an important point since the localization length 
cannot be defined through the exponential spatial 
decay of the electronic wave function in the finite 
aperiodic or disordered superlattices. This fact is 
a consequence of the strong spatial fluctuations of 
electronic wave function amplitude and therefore a 
more adequate description is based on the spatial 
extension of the wave functions given by the IPR 
parameter [4UI I5H 155] . 



IPR„ = / dx\ip n (x)\ 4 
h 



(10) 



A combination of equations ^ , ^ and ([9| leads 
to the relation between the generalized dimension 
of the electronic wave function and the scaling ex- 
ponent, T n (q), for the q-ih moment of the probabil- 
ity measure |43j . namely 



T n (q)=D n (q)(l-q). 



(11) 



Deviation of the scaling exponent from the linear 
function of q signals the multifractality of the elec- 
tronic state. 

In the limit of weak electric field, the interband 
transitions can be neglected and the influence of the 
electric field on the electronic states of Fibonacci 
superlattice can be considered by including an ad- 
ditional perturbation in the form 



W = -eFx . 



(12) 



where e is the elementary charge, and F is the ex- 
ternal homogeneous electric field applied along the 
growth axis of the system. 



In fact, the weak electric field is treated non- 
perturbatively as a parameter which produces the 
motion of the energy levels as it is shown in Fig- 
ure [2] For the considered superlattice the level re- 
pulsion leads to the formation of anticrossings in 
the energy spectrum, which is a simple consequence 
of the dimensionality of the system and the Dirich- 
let boundary conditions applied to the Schrodinger 
equation. 

Introducing the perturbation to the Schrodinger 
equation ([2| finally leads to the equation 



2m* dx 2 



tp(x) + [V(x) - eFx]rp(x) = EiP(x), 

(13) 

which allows us to investigate the effect of weak 
electric field on the global (energy spectrum, DOS) 
as well as local properties (LDOS, spatial fluctua- 
tions of the wave function amplitude) of the finite 
Fibonacci superlattice. 



3 Numerical results and dis- 
cussion 

Using the model of the finite Fibonacci superlattice 
described in Sec. [2] we have investigated its global 
as well as local electronic properties in the limit of 
weak electric field. All numerical values of the phys- 
ical quantities which are considered here are given 
in the atomic units, i.e., h = |e| = too = 1< The 
parameters of the model potential defined in equa- 
tion are taken to be V a = 0.27 eV = 0.01 a.u., 
c = 1.5 nm = 28.3 a.u., and p = 10, which means 
that the system consists of Alo.3Gao.7As barriers 
having width of 3 nm and separated by GaAs lay- 
ers of 3 nm or 6 nm width (the latter in case of 
two consecutive zeros in the Fibonacci binary se- 
quence). The number of barriers is N = 100 and 
the constant effective mass approximation is used 
with the value of m* equal 0.067too- 

One of the characteristic features of the obtained 
energy spectrum is the presence of minibands with 
fractal structure (Fig. [3]) . When the electric field is 
increased these minibands become broader, with a 
nonuniform distribution of the energy levels occur- 
ring during the entire process. It is the reason why 
darker and lighter areas corresponding to dense and 
rare subsets of the levels are present within the 
bands in Figure [3] showing the influence of the elec- 
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Figure 3: (Colour online) Density of states (DOS) 
as a function of energy E and electric field F. The 
results for the Fibonacci superlattice with N = 100. 
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trie field on the DOS of the Fibonacci superlattice. 
On the other hand, the general structure of the 
LDOS is not altered and the effect of the electric 
field is restricted to the change of the slope along 
the horizontal axis (see Figure|4]presenting how the 
electric field modifies LDOS, which leads to the pre- 
viously mentioned broadening of the minibands). 

The presence of the anticrossings between the en- 
ergy levels is the direct cause of the nonuniform 
structure of the energy spectrum as well as the den- 
sity of states function 35J. In the region of anti- 
crossings, the electronic wave functions change the 
degree of localization measured by the IPR param- 
eter. Figure [5] shows the IPR parameter for elec- 
tronic eigenstates of the Fibonacci superlattice as 
a function of electric field. The lowest degree of 
localization is observed in the absence of the elec- 
tric field, but in contrast to finite periodic systems 
the localization degree is not simply increasing with 
the electric field. In case of the finite Fibonacci 
superlattice the observed behaviour is much more 
complex. Instead of the successive states with very 
similar IPR dependence on the electric field, we can 
notice that for a chosen value of the electric field the 
highly localized states are separated by the states 
with much weaker localization. 

Therefore, the electric field affects the degree of 
localization, which in turn results in the change of 
the spatial fluctuations of the wave function. In 
order to explain the relation between the increas- 
ing electric field and the modification of the spa- 



Figure 4: (Colour online) Local density of states 
LDOS for the Fibonacci superlattice with N = 100 
as a function of energy E and position x. (a) F = 0, 
(b) F = 10" 7 [a.u.] and (c) F = 3 x 10~ 7 [a.u.]. 
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Figure 5: (Colour online) Inverse participation ra- 
tio IPR calculated for n = 0, . . . , 136 (i.e. for all 
eigenstates ip n which are bound for F = 0) in case 
of the Fibonacci superlattice with N = 100. 
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Figure 6: (Colour online) Differences (E n+ i — E n ) 
between the energies of the four lowest subsequent 
levels (upper panel) and the corresponding values of 
the inverse participation ratio IPR (lower panel) in 
the case of N = 100 Fibonacci superlattice. Letters 
(a)-(f) denote the positions of anticrossings. 



tial fluctuations of the wave functions we have per- 
formed a detailed analysis of the lowest energy lev- 
els, i.e. the ground-state level and the first two 
excited state levels. It should be noted, that each 
level is changed not only by the external electric 
field, but also by the coupling to the nearest en- 
ergy levels. As a result, a nonlinear dependence of 
the energy levels on the electric field is observed. 

The distances between the energy levels n and 
n + 1, calculated for n = 0,1,2, are presented in 
Figure [6] The neighbouring states tend to change 
their positions in such a way that the distances be- 
tween them increase or decrease alternately. This 
kind of analysis may be also generalized for higher 
energy levels, however, in this case a large number 
of anticrossings is observed which makes the anal- 
ysis not particularly suitable for the purpose of the 
clear explanation of the phenomenon. For this rea- 
son the further discussion focuses on the regions of 
anticrossings labelled by (a)-(f) in Figure [6j 

Figure [6] presents also the values of IPR param- 
eter calculated for the corresponding states in the 
same range of the electric field. For the values of 
the electric field at which the anticrossings are ob- 
served, i.e. values of (E n +i — E n ) have minima, 
also the values of IPR parameter for n and n + 1 
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Figure 7: (Colour online) Upper panel: Multifrac- 
tal spectrum f(a n ) calculated for the ground level 
state n = at different values of the electric field 
F with a m in and a max indicated for F = 0. Lower 
panel: The difference between the maximal and 
minimal singularity strengths, a max — a m { a , for the 
first four lowest energy levels, n — 0,1,2,3, with 
letters (a)-(f) denoting the positions of anticross- 
ings (see also Figure [fif. 



have local minima. It means that the electronic 
wave functions matching those levels become less 
localized and occupy a larger region of space due 
the coupling between the closely lying levels. 

The further study of the properties of electronic 
wave functions corresponding to the chosen states 
is based on the multifractal analysis. The values of 
the singularity strength a n and the corresponding 
singularity spectrum f(a n ), parametrized by the 
electric field, defined in equations Q and ([8]), are 
calculated using the standard 'box-counting' pro- 
cedure repeated for the values of the electric field 
changing from up to 10 _8 a.u. 

The f{a. n ) spectrum for the Fibonacci superlat- 
tice is found to be strongly asymmetric when no 
electric field is applied, as it is shown for n = in 
Figure [7J However, while f(a n ) is broadened when 
the electric field is present, the minimal value of 
the singularity strength, a m i n , remains almost con- 
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stant and close to zero, which is its minimal pos- 
sible value. This fact may be connected with the 
nature of the discussed eigenfuctions, which have 
the localized peaks almost not affected by the in- 
creasing electric field, as it can be seen in Fig. |4j 
Similar behaviour is observed also for higher eigen- 
states (see the lower panel of the Figure [7]), and 
resembles freezing of the a m i n value occurring in 
disordered systems [56] . 

As a result, the value (a max — c^min) is not an ob- 
jective measure of the electronic wave function lo- 
calization [57] in the case of this type of Fibonacci 
sequence-based systems under influence of the elec- 
tric field, which is clearly visible from the compar- 
ison with the values of the IPR parameter calcu- 
lated for the same states (Fig. [6]). Moreover, the 
first derivative of (a max — o> min ) is not continuous 
at the values of the electric field where the anti- 
crossings between the energy levels are observed. 
In light of these results, we rather think that for 
this kind of systems, the value (a max — QWn) can 
be used to detect the local minima of the IPR pa- 
rameter or equivalently to detect the anticrossings 
in the energy spectrum. 

The multifractal character of the electronic wave 
function in the finite Fibonacci superlattice and 
the influence of the weak electric field becomes 
much more visible in Figure [5J where the q de- 
pendence of the scaling exponent, T n (q), for a few 
lowest electronic states is presented. In all cases 
the scaling exponent is a convex function, mono- 
tonically increasing with q. The electric field mod- 
ifies the slopes of the considered scaling exponents: 
to a larger extent for negative q, whereas for posi- 
tive values of q the changes of the slope are much 
smaller. This type of changes of the slopes is related 
to marginal values of the singularity strength [2], 
amin (being almost constant as a function of the 
electric field) and a max (increasing notably in the 
electric field). 

Changes of D n (q) in the electric field are also 
noteworthy, as it is shown in Figure [9j In general, 
D n (q) decreases nonlinearly with q and reaches a 
constant for q — > — oo, however the value strongly 
depends on the state number on the other hand. 
For q — » +oo, D n (q) tends to a constant value, but 
the dependence of this value on the electronic state 
is rather weak. The nonlinear decrease of D„ (q) 
with increasing q is a presage of the wave function 
multifractality [43 . 




Figure 8: (Colour online) The scaling exponent 
r n (q) for the electric field F = (solid line) and 
F = 2 x I0~ 8 (dashed line) and for states n = 
0,1,2,3. 



In the following discussion we present results ob- 
tained for I < q < 3 which is the range where the 
fluctuations are mostly pronounced. Figure|9]shows 
the generalized dimension calculated in the vicinity 
of the anticrossings marked in Figure [6j D n (q) de- 
creases with q, but it has maxima when analyzed 
as a function of the electric field. The maximal 
values of D n (q) appear at the same values of the 
electric field as the anticrossings and are accompa- 
nied by the minima of IPR parameter. It leads to 
the conclusion that the maximum of the generalized 
dimension D n (q) as a function of the electric field 
corresponds to the decrease of the wave function 
localization degree. Moreover, if we analyze the 
surfaces of D n (q) for both states involved in the 
anticrossing and plot them as functions of q and 
the electric field F, it turns out that they cross and 
this crossing does not take place at the same value 
of the electric field for all values of q. The curve 
defined by the crossings that joins points for which 
the generalized dimensions are the same for the two 
discussed states, corresponds to the equality of the 
generalized moments fi n k of the probability densi- 
ties of the electronic states which are equivalent to 
the intensity fluctuations [55] . 
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4 Concluding Remarks 
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Figure 9: (Colour online) Generalized dimension 
D n (q) calculated in the vicinities of the anticross- 
ings (a-f) shown in Figure [6] 



In the present paper, we have investigated the influ- 
ence of the weak homogeneous electric field on the 
one-electron states in the finite superlattice gener- 
ated according to the Fibonacci sequence of two 
types of semiconductor layers. For clarity of pre- 
sented analysis only a few lowest-energy eigenstates 
have been considered in details, but conclusions can 
be generalized for higher-energy eigenstates. 

We have shown that the nonlinear character of 
the energy levels dynamics results from the large 
number of anticrossings. Therefore the energy 
spectrum of the Fibonacci superlattice and the den- 
sity of states are nonuniform in the limit of weak 
electric field. In the vicinity of the anticrossings, 
the inverse participation ratio of the electronic wave 
functions for individual electronic eigenstates pos- 
sesses minima. These minima of the inverse partic- 
ipation ratio are related to the maxima of general- 
ized dimension calculated as a function of electric 
field, that are observed for < q < 3 in the vicin- 
ity of the anticrossings. Moreover, we have shown 
that the change of the spatial extent of the electron 
wave function in the vicinity of the anticrossings 
is preceded by the correlation between the spatial 
fluctuations of wave functions corresponding to the 
electronic eigenstates participating in the anticross- 
ings. 

The relation between the positions of anticross- 
ings and the properties of the singularity spec- 
trum has been revealed, and the difference between 
the maximal and minimal values of the singularity 
strength has been found not to be a proper measure 
of the spatial extents of electronic wave functions 
in the presented case. Quite surprisingly, we have 
found the discontinuity of the first derivative of the 
difference between marginal values of the singular- 
ity strength for the values of the electric field where 
the local minima of inverse participation ratio ex- 
ist. We have also shown that a strong asymme- 
try of the singularity spectrum is preserved for all 
considered values of the electric field, although the 
electric field modifies the shape of the spectrum. 
This analysis of the relations between the singular- 
ity spectrum and the scaling exponents, together 
with the calculated generalized dimensions, have 
allowed us to show that the multifractal character 
of the electronic wave functions in this type of Fi- 
bonacci sequence-based systems is not destroyed by 
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the weak electric field. 

The results presented above correspond mainly 
to one arbitrarily chosen length of the Fibonacci 
sequence used as a basis for the model potential. 
Our further numerical studies performed for several 
different sizes of the systems show that the same 
characteristic features are observed, and thus we 
have decided to choose a typical example allowing 
the detailed analysis. 

We hope that the presented results can be useful 
in a deeper understanding of the nonlinear proper- 
ties of the energy spectrum in the aperiodic pho- 
tonic, phononic or semiconductor superlattices un- 
der influence of an appropriate external perturba- 
tion and motivate the experimental verification of 
the results. 
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